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ABSTRACT 



Context. In our SCUBA survey of Perseus, we find that the fraction of protostellar cores increases towards higher masses 
and the most massive cores are all protostellar. 

Aims. In this paper we consider the possible explanations of this apparent mass dependence in the evolutionary status 
of these cores. We investigate the implications for protostellar evolution and the mapping of the embedded core mass 
function (CMF) onto the stellar IMF. 

Methods. We consider the following potential origins of the observed behaviour: dust temperature; selection effects in 
the submillimetre and in the mid-infrared observations used for pre/protostellar classification; confusion and multi- 
plicity; transient cores; and varying evolutionary timescales. We develop Core Mass Evolution Diagrams (CMEDs) to 
investigate how the mass evolution of individual cores maps onto the observed CMF. 

Results. Wc find that two physical mechanisms - short timescales for the evolution of massive cores, and continuing 
accumulation of mass onto protostellar cores - best explain the relative excess of protostars in high mass cores and the 
rarity of massive starless cores. In addition, we show that confusion both increases the likelihood that a protostar is 
identified within a core, and increases mass assigned to a core. Selection effects and/or transient cores also contribute 
to an excess of starless cores at low masses. 

Conclusions. The observed pre/protostellar mass distributions are consistent with faster evolution and a shorter life- 
time for higher-mass prestellar cores. We rule out longer timescales for higher-mass prestellar cores. The differences 
in the prestellar and protostellar mass distributions imply that the prestellar CMF (and possibly the combined 
pre+protostellar CMF) should be steeper than the IMF. A steeper prestellar CMF can be reconciled with the observed 
similarity of the CMF and the IMF in some regions if a second opposing effect is present, such as the fragmentation of 
massive cores into multiple systems. 

Key words. Submillimeter;Stars: formation;ISM: clouds 
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> 1. Introduction 



The stellar Initial Mass Function (IMF) may largely be 
determined in the prestellar phase of the evolution of 
dense cores, as the close agreement between the form of 
the prestellar core mass function (CM F ) and the IMF 
entic in gly suggests dMotte et al. _ 19981; Testi k Sargent! 



ried out a submm survey with SCUBA (lHatchell et al 
20051 l2007bl lal hereafter Papers I— III, see also Kirk et al 
20061 and the similar 1300/im survey with Bolocam by 



19981: 



igl 
J 



Johnstone et alJ 120001 120011: iMotte et alJ l200lF 



Johnstone et all 120061: iReid fc Wilson! 120061: iKirk et all 
20061: iNutter k Ward-Thompson! 120071 ; I Andre et al.l l2007t 
Alves et alJl2007HGoodwin et al.ll2007h . However, feedback, 
competitive accretion and varying evolutionary timescales 
may all have a role to play in the evolution of the cores. 
Clearly, it is important to study how the mass in dense 
starless cores relates to the mass in cores forming proto- 
stars (protostellar cores), and therefore what determines 
the mass available to each protostar as it grows in mass. 

There are now a number of molecular clouds for which 
the dense cores have been catalogued (from millimetre 
(mm) and submillimetre (submm) surveys) and classified 
(using the mid-infrared detections from the Spitzer Space 
Telescope). Perseus is one such cloud, for which we car- 
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Enoch et alJ I2006D . We have identified over 100 submm 
cores, each of which has been classified as starless or proto- 
stellar using Spitzer (|Paper UUorgensen et alj|200dl2007l) 
and molecular outflows (|Paper Ilfl ). 

An intriguing result in Perseus is that the highest mass 
cores are all protostellar (Fig. [J), and as one moves towards 
lower masses, an increasing fraction of cores appear to be 
starless. This leads to the obvious questions: why are there 
so few massive starless cores and why are the two distri- 
butions different? In this paper we address these questions. 
We start by first reviewing the results of the previous work 
on this cloud (Sect. (2), and consider the effects of the anal- 
ysis of the submm data, detection limits and multiplicity 
on the distributions. Then in Sect. 2] we discuss some evo- 
lutionary scenarios which might produce these effects and 
introduce core mass evolution diagrams to illustrate how 
core evolution maps onto the observed CMF. In Section [S] 
we explore the factors which affect the relationship between 
the CMF and the IMF before presenting the summary and 
conclusions in Section [5] 
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2. Results 

The 103 su bmm cores discussed h ere were identified using 
Clumpfind ([Williams et alj 11994 ) on the SCUBA 850/xm 
map of Perseus spatially prefiltered by removal of a 2' 
Gaussian smoothed background (primaril y to remo ve arte- 
facts from the image reconstruction; see iPaper III for de- 
tails). The identified cores have masses in the range 0.5- 
50 Mq, sizes ~ 0.1 pc (1' at the assumed 320 pc distance 
of Pe rseus), an d mean H2 densities of > 3 x 10 5 cm~ 3 . 

In lPaperTfl , we classified the cores using mid-IR detec- 
tions from Spitzer to identify protostars. There are slight 
discrepanc ies between our classifi cation and the embedded 
YSO list of lJprgensen et all (|2007f ). due to the different clas- 
sification methods used (in our case, IRAC sources within 
12" or the C lumpfind source radius of the su bmm peak 
(see IPaper Hf ). whereas Ijprgensen et al.l (|2007f ) looked for 
MIPS 24/im sources within 15" of the pe ak). In total, 13 
of our sources are classified differently by I Jorgensen et all 
(2007), but only one of our massive protostellar sources 
(M > 1 Mq) is reclassified as starless. The majority of 
the massive protostellar sources are well-known (L1448C, 
L1448 NW, NGC1333 IRAS 4A, HH211). Therefore the 
alternative classification does not significantly change the 
distributions. We found that the ratio of protostellar cores 
to starless cores increased with mass, with ultimately no 
starless cores at all at the highest masses (above 12 Mq, 
IPaper J . Fig. 3). This result is shown in a modified form 
in the top panel of Figure [ITSection EO)) . 

In IPaper III] we used an alternative classification of the 
nature of the cores based on the presence or absence of 
molecular outflows for a subsample of 51 sources, and found 
the same result. 

In fact these effects (the most massive cores are 
protostellar, with more starless cores at lower masses) 
have now been in dependently confirmed in Perseus by 
lEnoch et all f2 008) using masses based on Bolocam 1.3/im 
data ( Enoch et al.l 120 06) . and the alternative Spitzer 



identification (iJprgensen et all 120061 iRebull et alj [200 



Uorgensen et al.H2007HLai et al J 12008) . (Enoch et all (12006 



also finds the same holds for Serpens, though intriguingly 
not Ophiuchus, which shows more similar mass distribu- 
tions for starless and protostellar cores. The massive star- 
forming region Cygnus X shows a more ext reme version 
of th is effect, with no massive starless cores (jMotte et alj 
12001 . 



3. Assumptions and selection effects 

Below we consider possible selection effects and assump- 
tions which could influence the our results and the core 
mass distributions. 



3.1. Detection and Identification 

The submm continuum detections are flux-limited. The ob- 
servations are not sensitive to cores below a column density 
threshold of 8 x 10 22 cmT 2 {A v ~ 90) (the 5ct flux detection 
limit at 850/im). For a point source, this corresponds to a 
mass limit of 0.5 M P) (assum ing 10 K and opacity k 850 = 
0.012 cm 2 g _1 , see IPaper Ilf ). but more massive sources can 
go undetected if the flux is more extended. Below this, 
our core detection is incomplete, and deeper studies are 
needed to determine the form of the mass distribution 



(eg. the Orion metastudy by iNutter fc Ward-Thompsonl 
2007). The sensitivity and high anticipated source counts 
of the next generation of submillime tre surveys, such as the 
JCM T Gould's Belt Legacy Survey (jWard-Thompson et al] 
2007 ) and Herschel surv ey of nearby star-forming regions 
( Andre fc Saracenoll2005l ) will lay to rest incompleteness is- 
sues at the 0.1 Mq level in nearby clouds. But currently, 
with a relatively shallow column density detection limit, for 
many cores we must be only detecting the 'tip of the ice- 
berg' with much of the mass lying at lower, undetectable 
column densities. In addition, objects over 2' (0.2 pc) in 
size are significantly affected by the spatial filtering of the 
data needed to removed artefacts on large size scales. The 
impact of the filtering depends on the degree of central con- 
centration of a source with the less condensed sources being 
the most strongly suppressed. 

For both these reasons (high column density detection 
threshold and filtering of extended cores) we are better at 
detecting cores with a high degree of central concentration, 
so the observations may only be tracing the later sta ges of 
prestellar evolution (eg. iTassis fc Mouschov ias 2004). For 
cores smaller than ~ 2' in size, higher mass cores are easier 
to detect. A priori we might expect a selection effect which 
favours the detection of starless (and protostellar) cores 
with higher masses. So the lack of high-mass starless cores 
is not a consequence of our inability to detect them if they 
exist. The submm detection limits do not explain the lack 
of high-mass starless cores. 

3.2. Temperature 



In lPaperTH the masses of the cores were calculated from the 
850/im emission assuming a constant dust temperature of 
10 K, which was the highest estimate for Td consistent with 
all the spectral energy distributions (SEDs). However, pro- 
tostellar sources tend to have higher temperatures and so 
produce relatively more 850/im emission for the same mass 
of circumstellar material. Can the assumed dust temper- 
ature reconcile the masses of the starless and protostellar 
cores? 

In order to test this, we recalculate the masses for the 
sources assuming different dust temperatures. Recent stud- 
ies of NH3 found a small difference between gas temper- 
atures of 11 and 12 K for starless a nd protostellar cores, 
respectively (jRosolowskv et al1l2007l ). However, differences 
in the dust temperature between starless and protostel- 
lar cores are evi dent from continuum observations of th e 
SEDs of sources (|Hatchell et al.ll2007bt lEnoch et al.ll2008l ). 
The apparent difference between the gas and dust temper- 
atures, averaged along the line of sight, is presumably due 
to the different sensitivity of molecular line and dust emis- 
sion to the warmer inner regions of the sources. Modelling 
of prestellar and protostellar cores has found equivalent 
isothermal dust temperatures of Td us t = 11 K (starless), 
14 K (Class protostars) and 16 K (Class I protostar s) 
(|Evans et al.ll200lh IShirley et al 1120021 1 Young et aHl2003t ). 

We follow lEnoch et all (|2006f ) in assuming 10 K for star- 
less cores and 15 K for protostars, in order to ensure that 
we account fully for dust temperature differences. The mass 
distribution for our sources assuming 10 K and 15 K tem- 
peratures for the starless and protostellar sources respec- 
tively is sho wn in the upper panel of Figure [TJ Although 
reduced from lPaper III , due to the lower masses inferred for 
the protostellar sources, there is still a significant differ- 
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Fig. 1. Top: Histogram of masses of submm cores, assum- 
ing 10 K for starless cores (blue line) and 15 K for proto- 
stars (red hatched). The SCUBA completeness limits for 
point sources at each temperature are marked at the top 
of the plot. Bottom: As above with protostars separated 
into Class (red hatched) and Class I (orange filled). 



ence in the mass distribution of the starless and protostel- 
lar cores, with more protostellar sources at higher masses 
and more starless cores at lower masses. To quantify this, in 
the whole population 46% of sources are starless, but above 
3 M Q this falls to 31%, 20% above 10 Mq, and zero above 
15 Mq. A K-S indicates that the probability that these two 
samples are drawn from the same mass distribution is 4%, 
compared to 0.001% wh en a temp erature of 10K is assumed 
for all the sources as in iPaper l3 . Only considering sources 
above 1 Mq, the K-S test gives a probability of 8% that 
the distributions are the same. 
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Fig. 2. Average SED for all 34 Class sources overplot- 
ted with modified blackbodies at 10, 15 and 20 K. Flux 
points are the mean of all detections from Class sources 
at each wavelength, from (left to right) 2MAS S, Spitzer 
IRAC, Michelle, IRAS, SCUBA and Bolocam (see lPaper III , 
for details). The 70 and 100 micron points are strictly up- 
per limits because of potential confusion within the IRAS 
beam, but still constrain the maximum temperature to less 
than 20 K. The 2MASS datapoints are marked as upper 
limits as the high fluxes are likely due to confusing sources. 

The lower panel of Figure [T] shows the core mass distri- 
bution but with the protostars separated into Class and 
Class I e volutionary stages based on the classification in 
IPaper iJ . As expected, the more evolved Class I cores have 
lower masses than the Class sources. From this figure we 
see that the excess of protostars is mainly due to the ex- 
cess of Class cores, and therefore it is the temperature 
of the Class population that is critical to the difference 
between the prestellar and protostellar mass distributions. 
Since the sense in which the mass distributions differ is that 
the Class sources are more massive, this difference could 
be reduced if the Class sources had even higher temper- 
atures. Raising the assumed temperature to 20 K for the 
Class sources would produce a distribution more com- 
patible with that of starless cores (15% K-S probability). 
However, a temperature as high as 20 K (which would be 
larger than the temperature of the Class I sources) is ruled 
out by the individual SEDs, in cluding those of the most 
massive Class cores ([Paper III ). The average SED for the 
Class cores is shown in Fig. [2] As the model SEDs show, 
the observations are compatible with a temperature of 15 K 
but not 20 K. Therefore even taking temperature into ac- 
count, there is still a significant difference between the two 
mass distributions. 

A consequence of assuming a temperature of 15 K rather 
than 10 K is that a given 850/im flux corresponds to a fac- 
tor 2 lower column density. As Clumpfind integrates mass 
down to a fixed flux limit, the mass of higher tempera- 
ture gas is measured down to a lower column density than 
lower temperature gas. If temperatures at the outer edges 



4 



Hatchell et al.: SCUBA Perseus survey - IV. Mass dependent evolution 



of protostars are significantly higher than those for starless 
cores, this creates a bias towards higher masses for proto- 
stars as surrounding material at lower column densities is 
included than for starless cores. To quantify this effect re- 
quires an understanding the temperature structure of the 
cores. The 15 K and 10 K values are the mass-weighted 
average in the beam, with the higher value for protostars 
due to the central heating. At large radii, the temperatures 
of protostars and starless cores are more similar. Models 
suggest that temperatures fall to below about 12 K within 
10,000 AU (30") of the central protostar and both starless 
and protostellar cores are heated to 12- 14 K at their outer 
edges by the inter stellar radiation field (jEvans et al.l f2001: 
IShirlev et al.ll2002T ). suggesting that this possible bias is not 
likely to produce a significant systematic effect. 

In addition to the different mass distributions of the 
protostellar and starless cores, a slightly different question 
is whether there is a significant excess of protostars asso- 
ciated with higher mass cores. Above 10 Mq, only 2 out 
of 10 sources are starless. Applying the binomial distribu- 
tion, the probability of getting 2 or fewer starless cores by 
chance is 9.3%. Although not entirely impossible that this is 
a random fluctuation (and we are limited here by the small 
number of high mass objects in one cloud), the overpopu- 
lation of protostars at high masses is likely to be real. This 
is further supported b y the excess of hig h mass protostars 
also found in Serpens (jEnoch et al.ll2008 ). 



3.3. Source Confusion 

The absence of massive starless cores could result from 
source confusion, if our analysis erroneously fails to sep- 
arate multiple lower mass cores, one of which happens to 
contain a protostar. This could occur due to insufficient 
angular resolution, or overlapping cores along the line of 
sight. Both of these are potentially more significant issues 
in regions where cores and protostars are clustered, such as 
NGC1333 and IC348. 

The 10 most massive protostellar cores in our sample 
are all in clusters or groups (Bl, IC348, L1448, NGC1333). 
This is consistent with the expectation that massive stars 
form in clusters. We might also expect to find the most 
massive prestellar cores in clusters. Yet only two of the 
four most massive prestellar cores are in NGC1333; the 
other two are the filaments to the SW of Bl and the SE of 
L1448, outside the main clusters of protostars. It is possi- 
ble that we fail to identify starless cores in clusters because 
they are confused with existing protostars. Many of the 
most massive cores contain multiple protostars: L1448 NW, 
NGC1333 IRAS 4A and 4B, NGC1333 IRA S 2A/B. The 
core Bl-b contains multiple submm peaks (|Hirano et al.l 
fl999h but only one object is identified as a protostar. 

In clusters, it is often unclear how to share the mass 
in the dense, fragmenting filaments between individual ob- 
jects. Multiplicity clearly affects our ability to identify, at 
the prestellar and protostellar stage, the amount of mass 
which ultimately will contribute to one star. As protostars 
produce strong peaks, and all clumpfinding software looks 
for peaks, some prestellar material may be being attributed 
to nearby protostars, increasing the mass of the protostellar 
cores at the expense of starless cores. Thus, confusion could 
explain why we identify the most massive protostellar cores 
in clusters. However, it does not explain why only the more 
massive starless cores are underabundant with respect to 



protostars, as we would expect at least as many lower-mass 
starless cores to be lost due to confusion. 



3.4. Stellar Content 

The difference between the mass distributions also depends 
on our ability to differentiate protostars from starless cores. 
For protostars, lumi nosity and outflow power decrease 
with decreasing mass (jBontemps et al.ll 199a [Fuller fe LadH 
I2002T ) . and therefore both the MIR and CO outflow emis- 
sion become weaker. The mid-IR and outflow detection lim- 
its could therefore explain why there are more apparently 
starless cores at lower masses. 

The Spitzer c2d survey (jEvans et al.ll2003T ) is limited in 
the lowest luminosity source it can detect. The Spitzer de- 
tections for the lowest-mass sources identified as protostars 
are all close to the IRAC detection limits and the addition 
of MlPS-only detections increases the num ber only by a few 
()Rebull et al.ll2OO7l ; lJ0rgensen et al.l2007|) . Spitzer c2d is es- 
timated to be 100% comple te to a luminosity of 0.3 Lq in 
Serpens (|Harvev et al.ll2007t ) though signifi cant fractions o f 
lower-luminosity sources are detected fe.g. lDunhamll2007l) , 
with 50% completeness at 0.01 Lq. We can assume these 
limits also apply to Perseus, as it is at a similar distance. 
So the population of protostars in Perseus with L < 0.3L Q 
remains incompletely sampled by Spitzer. For reference, a 
luminosity limit of 0.3 Lq implies completeness to sources 
ab ove 1 Mq final ste llar mass, according to the PMS tracks 
of ISiess et al.l (|1999D (assuming constant accretion). Some 
of the low-mass starless cores could contain undetected low- 
luminosity MIR sources and this will only be resolved by 
further deep searches in the MIR, such as with Herschel 
(lAndre fe Saracendl2005h . 

In lPaper IlJ we considered in detail the detection statis- 
tics for outflows and concluded that the outflow observa- 
tions also do not conclusively rule out the presence of low 
mass, low luminosity protostars in the apparently starless 
cores. Subsequent observations suggest that the majority 
of the cores are indeed starless, i.e. they do not contain 
outflows, when a deeper (by factor > 5) outflow search is 
carried out (Hatchell et al. in prep.). 

Nonetheless, although incompleteness might explain an 
excess of starless cores, it can only increase the fraction 
of starless cores as compared to protostars. In particular 
the failure of MIR and outflow observations to detect low- 
luminosity sources cannot explain why all the most massive 
cores are protostellar. 

3.5. Transient cores 

A population of transient cores, cores which are not gravi- 
tationally bound, might contribute to explaining the over- 
population of low-mass starless cores. A population of failed 
cores which will not form protostars or transient starless 
cores whic h incre ases at lower masses was suggested by 
Elmcgrecn ( 1997 ) and is seen in turbulence simulation s 
(jPadoan fe Nordlundl2002tlVazquez-Semadeni et al.ll2005l ). 
Cores whose chemistry su ggests transience are also ob- 
served (jMorata et al.ll2005l ). 

Estimates of the fraction of cloud mass in dense 
cores range from a f ew percent dFuller fe Mverd 119871 : 
Uohnstone et alj 12004 iMotte et al.l Il998f) to a few tens 
of percent ( Paper ll h typically with uncertainties of more 
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than a factor of two. Ultimately, the mass expected to end 
up in stars is also few percent of the total cloud mass 
(jLada fc L ada 2003), suggesting that in some cases there 
will be little mass in the core budget which does not ulti- 
mately go in to stars. 

It has been argued that cores that are centrally con- 
densed enough to be detected in the subm m must be form- 
ing s tars and therefore are not transient (jMotte fc Andrei 
1200 ID . With high mean densities above n# 2 ~ 3x 10 cm , 
corresponding to a thermal Jeans mass of 0.3 M Q ), our 
sources contain from one to tens of Jeans masses and thus 
satisfy the minimum condition to be gravitationally bound. 
On the scales probed by the N2H + line (critical density 
10 5 cm~ 3 ), lincwidths are roughly thermal and the majority 
of cores are virially supported if an external pressure term 
is included (jKirk et al.ll2007al ). It therefore seems unlikely 
that we are seeing a large population of unbound objects. 
However, it is possible that some cores could be destroyed 
by future energetic events, particularly cores in disruptive 
cluster regions, or that the currently binding external pres- 
sure decreases to the point where cores re-expand, particu- 
larly for lower-mass cores. 

4. Core mass evolution 

The observed prestellar and protostellar core mass func- 
tions are built up of mass measurements of many individ- 
ual sources. To fully understand them we must consider 
how individual objects evolve and how they appear to the 
observer at different times. By modelling the mass evolu- 
tion of cores, and reflecting this in the CMF, we have a 
tool to investigate how different factors in core evolution 
can lead to differences in the observed CMFs. 

One possibility is that the timescale for forming a star, 
i s f, may be dependent on core mass. As the number of 
prestellar or protostellar objects is the product of the pro- 
duction rate and the lifetime, longer lifetimes in a particu- 
lar phase imply a greater number of detections and a large 
representation on the CMF. 

A second factor is the relationship between the time it 
takes a core to form a star, and the time it takes for the core 
to grow in mass. The important parameter here is the ratio 
the timescale for a core of a given mass to form a protostar 
(t s f) to the timescale for the core mass to grow through 
the accretion of additional material from its surrounding 
cloud, (tgr). In the limit i s f Cf gr , a core can be considered 
as being fixed mass reservoir while it is forming stars. 

On the other hand, if i s f Z t gT , star formation will be 
taking place within a core as the core continues to grow in 
mass. These two scenarios have different implications for 
the relationship between the mass distributions of starless 
and protostellar cores. 

4.1. Core mass evolution and the CMF 

The detectable mass of starless cores and protostars 
changes as they evolve, so that each object moves across a 
range in observable masses during its lifetime, and the CMF 
is an average of sources at different stages in their evolution. 
Prestellar cores start as diffuse objects below the column 
density detection limits of submm surveys such as this one 
(see Sect. I3.1[) . and become more centrally condensed with 
time so that the detectable mass starts at the detection 
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threshold and rises throughout the prestella r phase. For ex - 
ample, in the seminal collapse calculation of lLarsonl ((T969), 
detectable central densities of 3 x 10 5 cm" 3 are reached af- 
ter 2 x 10 5 yr and the mass above the detection threshold 
continues to rise until a central core is formed at 5 x 10 5 yr. 
For protostellar cores, outflows and accretion remove mass, 
and the Class I sources have lower masses than Class 0s 

(Fig. ED. 

To represent this graphically we can plot Core Mass 
Evolution Diagrams (CMEDs, Figs. [3] and [4J which show 
mass vs. time for various models of core evolution, and the 
resulting core mass functions. For example, in Fig. we 
plot the evolution of a core with a simple linear evolution 
of mass with time, given by 

m(t) = m pcak (l - \t\/T), (1) 

where m pea k is the maximum mass reached by the core, t = 
is the time at which a protostar forms, and T = 0.5 Myr 
the timescale for core growth or mass loss. This repre- 
sents the accumulation of mass above the detection thresh- 
old during the prestellar phase, and loss of mass through 
accretion and outflow as a protostar. For this model, the 
time for core growth and the time for a star to form are 
equal: t s { = t gr . If we take a population of such cores and 
observe them at random times then we generate the CMD 
shown in the lower panel. The CMD is a power law with 
a positive index of +1, reflecting the longer times spent in 
the higher mass bins (which are equal width in log 10 (m)). 
Likewise, Fig. [3}d shows cores with evolution given by 

m{t) = m pcllk {\t\/T+l)- 1 ^ (2) 

with 7 = 1.35 and T = 10 4 yr. The resulting CMF is a 
power-law with negative index —1.35. This function, with 
longer timescales at lower masses, is the pathological case 
where time evolution alone reproduces a Salpeter-like CMF 
from a single population of identical sources. This mass 
evolution function becomes unrealistic at low masses be- 
cause the timescales tend towards infinite (note the ap- 
parent sharp cutoff at log lo (M/M ) = —0.85 is simply 
the mass reached at the time limits of our calculation). 
Nonetheless this, or a similar form for the mass evolution of 
cores, is interesting as it produces a CMF with a power-law 
slope at the high mass end similar to the slope of the IMF. 
Also, it implies a protostellar mass loss rate which peaks 
early in the protostellar stage and then dec reases with time 
in line with decreas ing outflow mass loss (jBontemps et al.l 
(|1996f ). but see also lPaper ill . 

In such a model, all stars would form inside high mass 
cores, with low mass stars produced in multiples or with low 
efficiency. The strongest evidence against such a scenario is 
the existence of low-mass low-luminosity Class sources, 
some of which have small outflows su pporting their youth 
(|Hatchell et al.ll2007at lDunham|[2007t ). A single peak core 
mass model is clearly extreme. However, a population of 
low-mass cores on the CMD in itself does not imply that 
protostars form in low mass cores. 

To generate more realistic diagrams we assume a pop- 
ulation of cores with a distribution of peak (maximum) 
masses m pea k following a three-part power-law (the peak 
core mass function or PCMF): 

dA^(m pcak )/(ilog(TOp Cak ) txTOp C ( " k _1) (3) 
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Fig. 3. Core mass evolution diagrams (CMEDs) for single evolutionary tracks. Each top panel shows on a linear-log plot 
the evolution of the core mass (as measured above a detection threshold) as a function of time. The bottom panels show 
the corresponding core mass functions assuming a population of identical cores sampled at random times, (a): Linear 
core mass evolution, (b): The pathological case where a population of identical cores results in a Salpeter power law 
CMF simply through mass evolution governed by Equation [2l 



with 

a = 0.3 0.2 < m poak < 0.8 

a = 1.3 0.8 < TOp 0ak < 5.0 
a = 2.3 5.0 < m P eak 



(-0.7 < log 10 m P oak < -0.1 
(-0.1 < log 10 TO pO ak < 0.7) 
(0.7 < log 10 m poak ). 



The boundaries between the power law segments are 
set a factor of 10 higher in m ass than those for the 
Galactic field IMF £K roupall20"02T) to take into account the 
star formation eff iciency ( Nutter fc Ward-Thomp son 2007; 
lAlves et al.l 120071 ) . The peak core mass function is signifi- 
cant because it gives the total mass available to each proto- 
star for accretion, which relates directly to the IMF through 
the SFE (see Sect. [5]). 

We then generate CMFs based on this distribution un- 
der various assumptions about the evolution of the core 
mass with time. From the core mass evolution we calculate 
the core mass function assuming that the ensemble of cores 
is observed at random points during their evolution. The 
observed population on the CMF in a given mass range is a 
sum over the population of cores of the time spent in that 
mass range. 

In the simplest case (Fig. |4ji), we assume that the mea- 
sured mass above the detection threshold increases at a con- 
stant rate in the prestellar phase, and decreases at a con- 
stant rate in the protostellar phase, following Equation!] 
The timescales for core growth and protostellar formation 
are equal, t gI = t s { = T . For this simple case the observed 
CMF is similar to the underlying peak mass function in 
that it shows a 3-part power law above 0.2 M Q . However, 
the power-law indices of 0.25,-0.4 and —1.35 are steeper 
than those of the underlying peak core mass function (0.7,- 
0.3,-1.3 in logm pea k) because of the contribution of evolv- 
ing cores detected at less than their maximum mass. This 
contribution increases towards lower masses, particularly 
where the underlying mass function has turned over, and 



also explains the population below log 10 m = —0.7 with 
index +1 which does not exist in the PCMF. 

Using CMEDs we can now investigate the effects of dif- 
ferent evolutionary scenarios on the observed CMF. 

4.2, Fast evolution of massive cores 

The lack of massive starless cores, and the decrease in the 
number of starless cores towards higher masses, could im- 
ply that more massive cores have relatively short prestellar 
lifetimes, and that this is why we see few of them. (The 
alternative, that higher-mass protostars have longer life- 
times, can not explain why so few massive starless cores 
are detected). Current best estimates for the mean lifetime 
for starless cores detectable in the submm are similar to 
the protostellar lifetime, 1.5-4 x 10 5 years (|Hatchell et al.l 
120051: lEnoch et al.l 120071 iKirk et all l2007bf) . but the free- 
fall timescales for mean densities of > 3 x 10 5 cm -3 are 
less than 7 x 10 4 years, so theoretically this observed phase 
could be very short. 

From the observations we can estimate the prestellar to 
protostellar ratio for cores at masses above 10M Q : there 
are four times as many protostars suggesting a prestellar 
lifetime 1/4 of that in the protostellar phase. For compar- 
ison, at 3 Mq the number of cores in the prestellar and 
protostellar phases are roughly equal. 

In Fig. QJ) we plot a CMED in which massive cores 
evolve fast in the prestellar phase, with mass evolving lin- 
early with time according to Equation [T] but with varying 
prestellar timescale T = 1/ i/wpcak/ 'Mq Myr. The proto- 
stellar timescale remains constant at T = 0.5 Myr. The re- 
sulting prestellar mass distribution is steeper than the core 
peak mass function (steeper than to -135 ) and steeper than 
the (unchanged) protostellar mass distribution. Therefore 
we can reproduce the excess of high-mass protostars by as- 
suming protostellar lifetimes are all equal but the prestellar 
lifetimes depend strongly on mass, with short lifetimes for 
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Fig. 4. Core mass evolution diagrams (CMEDs) for different evolutionary scenarios. Each top panel shows the evolution 
of core masses as measured above a detection threshold as a function of time for various peak masses m pea k, on a linear- 
log plot. The bottom panels show the corresponding time-averaged core mass functions for prestellar cores (blue, vertical 
hatching), protostellar cores (red, diagonal hatching), and all cores (thin solid line). The straight lines are marked with 
approximate power laws for sections of the CMF. (a) Equal pre/protostellar timescales for all cores regardless of mass, 
with a constant mass accumulation/removal rate ±m oc m pca k. (b) Faster evolution for more massive sources in the 
prestellar phase only, (c) Continuing accretion into the protostellar phase, so that the peak mass is reached some time 
after t = 0. 



high mass sources. Note that this model still forms stars at 
the point at which the core stops growing, t s { = t gT . 

Longer prestellar lifetimes for lower-mass cores 
are predicted by am bipolar diffusion models 
(jTassis fe Mouscho vias 2004), and our results are qualita- 
tively consistent with this. Starting from the assumption 
that higher mass cores evol ve from lowe r dens ity (hence 
have a higher Jeans mass), IClark et alj (|2007h find the 
opposite: that the free-fall timescale is longer for higher 
mass cores, and therefore the prestellar mass distribution 
should be flat ter than the I MF. O ur result differs from that 
presented by IClark etldl (|2007l ) because our cores evolve 
in mass during the prestellar and protostellar phase, so 
there are differences in how the observed CMF is derived. 
Nonetheless, longer timescales for higher mass clumps are 
inconsistent with our results. 

The above discussion assumes that cores evolve at fixed 
(gas+stellar) mass, ie. the timescales for star formation are 
much less than the timescales for acquisition of additional 
material by the cores, t s { <C t gr . In this case the difference 
in mass distributions argues that massive cores must them- 
selves form and then evolve to form stars more rapidly than 
low mass cores. 

4.3. Continuing accretion 

In a model where the core is not considered to be a fixed 
mass reservoir but continuing accumulation of mass is al- 
lowed during the protostellar phase, the envelope mass 
can continue to rise during the main protostellar accretion 
phase. This is the case where the star formation timescale is 
less than the timescale for core growth, f s f < t gI . The peak 
envelope mass can then be reached sometime during the 
protostellar phase rather than at the prestellar/protostellar 



transition. Such models require surrounding lower-density 
material to become gravitationally bound to the core and 
continue to flow onto the envelopes of existing protostars. 

A model of this is shown in the CMED in Fig. 2J;. In 
this case, the time at which the peak core mass is reached 
is shifted into the protostellar phase by an amount propor- 
tional to log 10 Wpoak — log io m irim where m min — O.2M 
is the minimum mass in the peak core mass function, 
so i gr > t s {. The resulting prestellar mass distribution is 
steeper than the peak mass function and steeper than the 
combined (pre+protostellar) CMF. The protostellar mass 
function is flatter, reflecting the reduced amount of time 
spent at low masses. Again, this model reproduces the 
main features of the observed pre/protostellar mass dis- 
tributions. Core masses which continue to rise into the pro- 
tostellar phase can explain why protostellar cores tend to 
have masses higher than those of starless cores. 

Core masses will continue to increase into the proto- 
stellar phase as long as the mass loss through accretion 
and outflows is outweighed by mass moving from low den- 
sity at large radii to within the column density detection 
threshold of the observations. Although the model in our 
CMED is clearly simplistic, there are several ways in which 
this kind of behaviour can arise. Recent models of core 
growth through ambipolar diffusion show significant in- 
ward velocities in the outer parts of the c ore which con- 
tinue after the central protostar has formed (jAdams fe Shul 
I2007D . Alterna tively, if large-scale turbul ent flows supply 
the cores (see iMac Low fe Kles scn 2004) then cores will 
continue to grow until the supply of material is exhausted. 
Additionally, in competitive accretion models, cores which 
have already accumulated higher masses form deeper grav- 
itational wells and competitively accrete gas from the sur- 
rounding cloud and lower mass cores, which might bias the 



8 



Hatchell et al.: SCUBA Perseus survey - IV. Mass dependent evolution 



grow th towards already massive protostars (jBonnell et al.l 
1200 ll but note that there are kinematic arguments against 
competitive accretion: [Walsh ct al. 2004-. jAndre et aO2007t 
J0rgensen et al.M2007l iKirk et al.M2007al though see also 



Avliffe et alJl2007ft 



4.4. Luminosity constraints on mass evolution 

The bolometric luminosity is plotted against core mass for 
the Perseus sources in Paper II, Fig. 4 (protostars) and 
Paper III, Fig. 2 (all sources). There is a general trend 
of increasing luminosity with increasing mass, but with a 
large scatter. There are sources with an order of magni- 
tude difference in core mass (the complete sample spans 
0.5-50 M ) that produce similar bolometric luminosities 
while there is an order of magnitude scatter in source lu- 
minosity over some ranges of core mass. In principle the 
observed luminosity-core mass distribution provides a con- 
straint on the evolution of the core mass and hence the 
CMF, but exploring this is complex. 

The mass evolution of cores affects the luminosity of the 
protostars forming in the cores through accretion of core 
material onto the central star-disk system. The accretion 
luminosity is given by L acc — GM*M acc / R acc where M* is 
the stellar mass, and M acc and R acc are the accretion rate 
onto the star/disk and the accretion radius respectively. 
However to predict the luminosity evolution of the sources 
requires an understanding of how the accretion rate onto 
the star/disk M acc and the rate of change of mass of the 
core m are related. This relationship is likely to be both 
complex and evolving. In addition to the mass loss due 
to accretion on to the central source, m encompasses all 
the other core mass loss, and growth, mechanisms. This 
includes, for example, the mass loss associated with the 
effect of winds and outflows (not necessarily self-generated) 
as well as the growth of the core due to the accretion of 
surrounding cloud material. So for example, in the case 
where a core continues to grow in mass once a protostar has 
formed, M acc and rh have opposite signs up until the point 
when the core ceases to grow in mass. Even in the simple 
case where M acc is a constant fraction of m, the radius 
into which the material falls, R acc , and hence L acc , is likely 
to evolve, as, for example, a circumstellar disk forms and 
grows, and core material of different angular momentum is 
accreted. 

Using various simplifying assumptions and initial con- 
ditions, this luminosity e volution has been modelled 
by a number of authors d Smith! l2000t Froebrich et alj 
2006; lMverdfl998t iSiess et alJll99fl ISaraceno et al.lll996h 
However these models have mostly concentrated on cores 
which contain a fixed reservoir of material and have not in- 
vestigated the evolution of cores as massive the most mas- 
sive cores observed in Perseus. It would be interesting to 
extend such models to encompass both the range of core 
masses seen in Perseus and possible alternative scenarios 
for the evolution of the core mass, but such models are 
beyond the scope of this paper. 



5. From CMF to IMF 

So far we have discussed how the observed core mass dis- 
tributions relate to the underlying peak core mass distribu- 
tion, and how differences between the mass distributions of 



the prestellar and protostellar cores might arise. The sub- 
sequent issue is how these cores go on to form stars with a 
range of masses consistent with the IMF. 

Observational selection effects and systematic er- 
rors (Section [51) , as w ell as small number statistics 
(Swift k Williams 12008!) . make current measurements of 
the CMF highly uncertain. Nonetheless, the apparent sim- 
ilarity of the power-law slope of the CMF and the IMF has 
now been noted in several regions. The simplest explana- 
tion for the match is that i) there is a one-to-one mapping 
between cores and stars which are formed, ii) there is a 
uniform star formation efficiency (SFE) for all the cores 
i.e. that the same fraction of the mass of each core is con- 
verted in to star, Hi) the measured mass of each core is 
representative of the total mass available to form a star, 
and iv) the evolutionary timescales are independent of the 
core mass. 

If any of these assumptions is not in fact true, then the 
mapping of the CMF to the IMF of the stars formed is 
clearly more complex. However we can consider breaking 
each assumption in turn and examine how the resulting 
IMF relates to the CMF as measured for prestellar cores 
and/or all submm cores. 

5.1. One-to-one mapping 

Multiplicity breaks the one-to-one mapping between cores 
and stars. This has b een c onsidered in detai l recen tly by 
iGoodwin et ail ([20071) and ISwift k Williams! (|2008fl . who 



demonstrate that if multiplicity is significant we would 
expect the CMF to be skewed to higher masses com- 
pared to the IMF, as massive cor es contribute to a fi- 
nal population of lower-mass stars. IGoodwin et alj (|2007f ) 
argue that the best match betw een the observ ed CMF 
dNutter k Ward-Thompson! 12007! ) and the IMF (jKroupal 
2002) occurs if all cores (including low-mass cores) form 
binaries or higher-order multiples. 

There is plenty of evidence for multiplicity in Perseus 
sources. Although multiple Spitz er MIPS sources ca n 
only be identified in 3 cores (jjprgensen et alj |2006[) . 
at least 10% of the SCUBA cores have been demon- 
strated to contain multiple sources when observed with 
high- resolution infrared, radio or mm/submm interferome- 
trv (lLav et al J 1 19951: lAnglada et aTl 120001: [Rodriguez et al l 
1 19971: iRodriguez k Reipurthlll998t I Wolf-Chase et al.ll2000t 
lO'Linger et al.ll2006tlHirano et al.lll999ft . and the majority 
of the submm sources have not been the targets of detailed 
studies of their fragmentation and multiplicity. 

A population of cores which will never form stars, tran- 
sient cores, would also break the one-to-one mapping be- 
tween observed cores and stars. Transient cores are more 
likely to appear at low mass, where the gravitational bind- 
ing is weaker, and skew the CMF to lower masses than the 
IMF. The presence of a population of low-mass, transient 
cores cannot be ruled out by our data. 

5.2. Star formation efficiency 

If the star formation efficiency (SFE) is not independent of 
core mass, variations in it can have a complex effect on the 
relationship between the CMF and IMF. Changes in the 
SFE can lead to an IMF which is either steeper or flatter 
than the CMF or affect different mass ranges differently, 
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depending on exactly ho w the SFE varies with core mass 
(jSwift fe Williams! l2008h . if indeed the SFE is an unique 
function of core mass at all. SFE variation therefore acts 
to skew or blur the relationship between the CMF and the 
IMF. 

Although there are estimates of the instantaneous 
star formation efficien cy in cores in Perseus of 10-15% 
(jJorgensen et al.ll2007f l. we have no information on the vari- 
ation in SFE between cores. Future constraints on the SFE 
may come from studying the mass transfer through accre- 
tion and outflows and probing how feedback processes affect 
star formation within a core. 

5.3. Representative mass measurements 

If the core masses measured are not representative of the 
total mass available for accretion over the lifetime of the 
core, then the relationship between the IMF and the CMF 
is again more complex. As discussed in Sec. 4.1, the ob- 
served mass of a core (or distribution of masses of a sample 
of cores) is a function of the peak mass of the cores and 
the time evolution of their mass. Depending on these two 
factors, different ranges of observed core masses may con- 
tribute to the same range of stellar masses in the IMF. In 
general cores will be observed at masses less than their peak 
mass, skewing the CMF to lower masses compared to the 
IMF. Only if cores spend a large fraction of their evolution 
with masses at, or near, their peak mass is the observed 
mass of core representative of the total mass available to 
form stars, leading to a one-to-one relationship between the 
IMF and CMF. 

5.4. Evolutionary timescales 

Longlasting evolutionary phases will be relatively overrep- 
resented in the count of cores on the CMF. Conversely, the 
most rapid phases are the least likely to be seen. Therefore, 
if cores of different peak masses evolve at different rates, 
the shape of the CMF is affected. This is in contrast to the 
IMF, which takes account of the short lifetimes of massive 
stars by means of stellar evolution theory. 

The evidence presented here for the absence of massive 
prestellar cores in the CMF can be explained if cores with 
high peak masses evolve more rapidly. Therefore, from the 
observed deficit of massive prestellar cores we would expect 
the IMF to be flatter (higher populations at higher masses) 
than the CMF. 

Putting all four factors together, what can we say about 
the relationship between the measured CMF and the IMF? 
It is clear that if high-mass prestellar cores are underrep- 
resented on the CMF, and cores map onto the IMF one- 
to-one with a constant core star formation efficiency, then 
we would expect the measured CMF to be steeper than the 
IMF at the high-mass end. At the low mass end, the CMF 
always contains an additional population of cores evolving 
to or from higher masses and should always be overpopu- 
lated relative to the IMF. To reconcile the observed similar- 
ity between the CMF and the IMF, it would be necessary 
to invoke two opposing effects, eg. faster evolution for more 
massive cores combined with multiplicity, which contrive 
to cancel out so that the observed CMF once again, and 
coincidentally, has the same power-law form as the IMF. 



6. Summary and conclusions 

The mass distributions of the prestellar and protostellar 
cores in Perseus are different, with an excess of protostel- 
lar cores at high masses and an excess of starless cores at 
low masses. There are a number of possible selection effects 
which could be influencing the observed distributions. Of 
these, only confusion in regions where cores and protostars 
are highly clustered can plausibly explain the observed dis- 
tributions. 

In clustered regions there is a difficulty in identifying 
the structures which will ultimately form a star. This may 
be because the gas which will go on to form future stars 
lies in dense, filamentary structure with no clear bound- 
aries segregating it into discrete, centrally peaked cores. 
Confusion in the observations then leads to both an un- 
dercounting of starless cores and also an overattribution 
of mass to protostars if they are confused with essentially 
prestellar material. 

Of the other assumptions and selection effects involved 
in deriving the mass distribution, differences in the temper- 
atures between classes of cores cannot explain the differ- 
ences between the mass distributions. At low masses, the 
distributions are affected by selection effects, particularly 
our increasing inability to identify low-luminosity sources 
as protostellar. Transient cores may also contribute to the 
excess of starless cores at low masses, but seem unlikely to 
be significantly influencing the relative underabundance of 
massive starless cores. 

If selection effects are not the most significant cause of 
the difference in mass distributions, then these differences 
must be direct consequence of the evolution of cores. By 
looking at the mass distributions of the prestellar and pro- 
tostellar cores separately, we can place constraints on core 
mass evolution. We find two evolutionary scenarios which 
explain the observed excess of protostellar cores with high 
masses. The simplest is that timescales for the formation of 
protostellar cores may vary with core mass, with more mas- 
sive cores forming protostars more quickly than lower mass 
cores. Alternatively (or additionally), the masses above our 
column density threshold may continue to rise well into the 
protostellar phase, rather than peaking at the point when 
a protostar is formed. The observed distributions rule out 
the possibility that higher mass cores evolve more slowly to 
form protostars than lower mass cores. 

Whatever the explanation for the observed differences 
in the mass distributions, they provide a clue to how cores 
evolve to form stars, and how the measured CMF relates to 
the stellar IMF. The relationship between the CMF and the 
IMF is in general complex and several assumptions about 
evolution, detectability and efficiencies have to be made to 
derive the IMF from the CMF. But clearly if we are under- 
counting the number of massive starless cores which form, 
either because they only last a short time or because of con- 
fusion with protostars, then the IMF should contain more 
systems at the high mass end of the distribution than the 
prestellar CMF. The IMF should also show relatively fewer 
sources at the low mass end compared to the CMF, because 
the CMF contains a contribution from cores observed while 
evolving to or from their maximum mass. On the other 
hand, we know that many cores will fragment into multipl e 
lower-mass systems on the IMF (see lGoodwin et afl l2007). 
In some clouds these effects may conspire to cancel each 
other, producing a CMF and an IMF with similar shapes. 
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Clouds such as Perseus where the distribution of prestel- 
lar and protostellar core masses are inconsistent with an 
apparently simple, direct mapping of the CMF to IMF pro- 
vide important insights into the possible evolution of cores 
as they form stars. Deeper extensive surveys of cores and 
protostars in other clouds are needed to help unravel the 
details of these possible evolutionary paths of the cores. 
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